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We calculate the amount of polarization-entanglement induced by two-photon interference at a 
lossless beam splitter. Entanglement and its witness are quantified respectively by concurrence 
and the Bell-CHSH parameter. In the presence of a Mandel dip, the interplay of two kinds of 
which-path information — temporal and polarization — gives rise to the existence of entangled 
polarization-states that cannot violate the Bell-CHSH inequality. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.50.Dv 
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I. INTRODUCTION 

Entanglement, the nonclassical correlations between 
spatially separated partides, is typically a signature of 
interactions in the past or emergence from a common 
source. However, it can also arise as the interference 
of identical partides pj. By postselecting experimental 
data based on the "click" of detectors @,U, photons scat- 
tered at a beam splitter have violated a Bell inequality, 
even if they originated from independent sources 0, Q . 
In reverse, triggered by an interferometric Bell-state mea- 
surement, entanglement has been swapped |fj to initially 
uncorrelated photons of different Bell pairs 0, @ • The 
observation of these nonclassical interference effccts is an 
important step on the road towards an optical approach 
of quantum information processing [Ï3, uM ■ 

Being furnished by interference, the ability of a beam 
splitter to entangle the polarizations of two independent 
photons depends on their indistinguishability 12]. One 
of the incident photons is horizontally polarized in state 
|H;?/>), the other vertically polarized in |V; 4>)- The pho- 
tons are partially distinguishable by their temporal de- 
grees of freedom captured in the kets and \d>). Besides 
temporal which-path information inherited from incident 
photons, a scattered two-photon state possibly holds po- 
larization which-path information. We make no assump- 
tions about the scattering amplitudes connecting polar- 
izations at the beam splitter, except that they constitute 
a unitary scattering matrix. Translated to a polarization- 
conserving beam splitter, this corresponds to incident 
photons in states |ct;ï/>) and \cr';(f>) where er, a' are ar- 
bitrary superpositions of H, V. Our analysis generalizes 
existing work on a polarization-conserving beam splitter 
where a = H and a' = V [lE^. 

The polarization-state p of a scattered photon pair is 
established from the scattering amplitudes of the beam 
splitter, the shape and timing of photonic wavepackets 
(\ip), \4>)) and the time-window of coincidence detection. 
If not erased by ultra-coincidence detection, an amount 



of temporal distinguishability of (1 — 1(^10) | 2 ) pertains 
corresponding to a mixed state p. We calculate both its 
concurrence and the Bell-CHSH parameter. The ability 
of the latter to witness entanglement can disappear in 
the presence of a Mandel dip. In terms of a polarization- 
conserving beam splitter, this corresponds to a deviation 
of <t, a' from a = H and a' = V. 



II. FORMULATION OF THE PROBLEM 

In a second-quantized notation, the incident two- 
photon state |H;i/>)l|V; </>)r takes the form 

l*in) = 4,L<r|0), (1) 

with field creation operators given by (see Fig. ^) 



*^ L = / dw<4(w)^*(w), $V,R .= / dub^uj)^*^). 

(2) 

(The subscripts R,L indicate the two sides of the beam 
splitter.) The operators ai{w) with i = H, V satisfy com- 
mutation rules 

[oíH, <»,■(</)] = 0. [oíH,o}(w')] - %<5(w - to'). (3) 

The same commutation rules hold for the operators bi(u>), 
with commutation among a and b. 

The outgoing operators Cí(uj), di(ui) are related to the 
incoming ones a,i(ui), bi(u>) by a 4 x 4 unitary scattering 
matrix S, decomposed in 2 x 2 reflection and transmission 
matrices r,t,t' ,r': 

c(u) \_frt'\( a(u) \ ( a H (w) \ 

d(w) )-{t r> ){ b{uj) ) ' a[UJ) ~ { ayiuj) ) ' 

and vectors í>(w), c(u>), d{uj) defined similarly. The scat- 
tering amplitudes are frequency-independent. 
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The outgoing state l^out) can be conveniently written in a matrix notation 



|*out> = du dw'^Wft/) 



ct( w ) 
dt( w ) 



rer in í' T rcr in r' T 



C t( W ') 



|0). 



(5) 



r 






c 









FIG. 1: Schematic drawing of generation and detection of 
polarization-entanglement at a beam splitter. The indepen- 
dent sources SL and SR each create a photon in modes {a} 
and {ò} cf. Eq. Q. The beam splitter with unitary 4x4 scat- 
tering matrix S couples the polarization of incoming modes 
to the polarization of outgoing modes {c} and {d}. Polariza- 
tions are mixed by Rl,Rr and absorbed by photodetectors 
DL,DR. A coincidence circuit C registers simultaneous detec- 
tion of photons. 



Here we used the unitarity of S and a- m — (a x + ia v )/2, 
with a x and a v Pauli matrices, corresponds to the polar- 
izations of the incoming photons cf. Eq. 0J. The matrix 
Cin has rank 1 reflecting the fact that polarizations are 
not entangled prior to scattering. Since we make no as- 
sumptions about the scattering amplitudes (apart from 
the unitarity of S), the choice of a m is without loss of 
generality (see Appendix[SJ. 

The joint probability per unit (time) 2 of absorbing a 
photon with polarization i at detector DL and a pho- 



ton with polarization j at detector DR at times t and t' 
respectively is given by 

Wij (t,t') <x (* «t|^ ) (t)£?fe ) (É')s5+ ) (í / )£;W(t)|* aut ), 



(6) 

where (t) and E^ (t) are the positive frequency field 
operators of polarization i at detectors DL and DR. The 
probability Cy (t) of a coincidence cvent within time- 
windows t around t is given by 



Ca(t) 



dt' / dt"wij{t',t"). (7) 



Experimentally, the time-window r has typically a lower 
bound determined by the random rise time of an 
avalanche of charge carriers in response to a photon ab- 
sorption event. 

The polarization-entanglement is detected by violation 
of the Bell-CHSH inequality This requires two local 
polarization mixers í?l and í?r. The Bell-CHSH param- 
eter £ is 

£ = \E(Rj j ,Rb 1 )+E(R 1j ,Rji) + E(Rj j ,R r ) — E(R l ,R k )\, 

(8) 

where E(Rl,Rr) is related to the correlators 
Cij{R h ,RR) by 



E 



Chh + Cw — Chv — Cvh 
Chh + C*vv + Chv + CVh 



(9) 



Substituting the correlators of Eq. Q into Eq. JSJ| , we 
see that 

E{R L ,R R ) = Ty P {RÍ<t z R l ) ® {R^Rr), (10) 

where a z is a Pauli matrix and p a 4 x 4 polarization 
density matrix with elements 



Pij,mn = jf((l+ |a| 2 )(7l)ij(7l)mn + i 1 ~ H^Mij (72) mn) ■ 

The parameter a is given by 



It-Ï dt' Jdiüf dcü'4>(u)i>* (w') 



*rw'v*( w - w ')*' 



J(ft-* dt> fduf ^'<K w )0*( w O ei(ü ' _ü '' )t ') (jl^l dt'fduj (L·'^(u))^*{üj')é^- U ') t '^ 



r 



(ii) 



(12) 



and 71,72 are 2x2 matrices related to the scattering amplitudes by 

71 = ra in r' T + t' V£t T ', 72 = ra in r' T - t> ' aJJ 1 '. (13) 
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The normalization factor J\í takes the form 

N = (1 + M 2 )Tr 7 í 7l + (1 - |a| 2 )Tr 72 t 7 2 . 



(14) 



The parameter 1 — \a\ 2 E (0, 1) represents the amount 
of temporal which-path information. Generally, the time- 
window r is much larger than the coherence times or 
temporal difference of the wavepackets. We may then 
take the limit r — > oo and a reduces to the overlap of 
wavepackets 



(15) 



In the opposite limit of ultra-coincidence detection where 
t — > 0, temporal which-path information is completely 

|2 



erased corresponding to \a\ = 1. 



III. ENTANGLEMENT OF FORMATION 

The entanglement of formation of the mixed state p is 
quantified by the concurrence C [l^ | given by 



(16) 



C = max ( 0, v Ai 



The Ai's are the eigenvalues of the matrix product pp, 
where p — (a y <g> a y )p*(a y ® o y ), in the order Ai > A 2 > 
A3 > A4. The concurrence ranges from (no entangle- 
ment) to 1 (maximal entanglement). For simplicity of 
notation it is convenient to define (xy)ij, m n = XiíUmn- 
The matrix p can be written as 



P = jf ((1 + H 2 )7i7i + i 1 ~ l"| 2 )7272 



(17) 



with 7 = a y j*a y . The product pp takes the simple form 



PP 



Tr7i7i 
N 2 



(1 + M 2 ) 2 7i 7 i - (1 - HT 7 272 ) 



2\2, 



(18) 



where we have used the multiplication rule xyvw 
(Tiy^v)xw and 



Tr7i7i 



-Tr 7 t 72 , Tr j\% = Tr 7 | 7l = 0. (19) 



The results for the tilde inner produets of Eq. (|19fl hold 
since the photons are not polarization-entangled prior to 
scattering (Det cr; n = 0). 



The non-Hermitian matrix pp has eigenvalue-eigenvector decomposition 

((1 



PP = 




\a\ 2 ) 2 s lSl 




(20) 



where we have defined orthonormal states si = (1/2)(1 + 
a z ) and S2 = {l/2){o'x+io'y)· The pseudo-inverse is easily 
seen to be 



It follows that 



\ — — (sxji - S2J2] ■ (21) 

(Tr 7 Í7i)* V ' 



2|«| 2 |Tr 7l t 7l 
M 



C 



(22) 



The trace that appears in the numerator of Eq. i|22[l is 
given by 



ITT7Í7J = 2^/DetXtXDet(l - X*X), (23) 
where we have defined a "hybrid" 2x2 matrix X as 

x = (ruu Jiív ^ (24) 

V r VH íyv ) 

The normalization factor M given by Eq. (|14(l can be 
expressed in terms of X using 



Tr7Í 7l = Tr X^X - 2PerAT t X, 



(25) 



Tr7| 72 = TrX+X - 2DetX t X. 



(26) 



( "Per" denotes the permanent of a matrix.) In the deriva- 
tion of Eqs. I|23I25I26I) we have made use of the unitarity 
of S. The concurrence becomes 



C 



2|a|VDetXtXDet(l - XVC) 



TrX^X - (1 + \a\ 2 )PcrXïX - (1 - \a\ 2 )DetXïX' 

(27) 

Entanglement depends on the amount of temporal in- 
distinguishability \a\ 2 and the Hermitian matrix 



yty - I l r n| 2 r H -tV 



(28) 



containing the states rn = (fHHj^vl·l) and t v = 
(íJj V ,í vv ) of a reflected and transmitted photon to the 
left of the beam splitter. The determinant of X'X mea- 
sures the size of the span of rn and ty as 



DetX f X 



|r H | 2 |t v | 2 



|r H -tyi 2 

|r H | 2 |tVI 2 



(29) 
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If i'n and ty are parallel (DetX^X = 0), a scattered pho- 
ton to the left of the beam splitter is in a definite state, 
giving rise to an unentangled two-photon state (C = 0). 
Similarly, 

Det(l - Xtx) = |t H | 2 |r v | 2 (l - (30) 

involves scattered states tn = (íhh,ívh) and r v — 
( r HVi r vv) t° the right of the beam splitter. The denom- 
inator of Eq. (|27() is the probability of finding a scattered 
state with one photon on either side of the beam splitter. 
It deviates from its classical value (X* X)hh+(^^)vv — 
2(X^X) Hïi (X J <X) vv by an amount -2|a| 2 |(XtX) H v| 2 
due to photon bunching. This reduction of coincidcnce 
count probability is the Mandel dip [ÏTj ■ It measures the 
indistinguishability of a reflected and transmitted photon 
as the product of temporal indistinguishability \a\ 2 and 
polarization indistinguishability |(XTJf )hv| 2 - 

IV. VIOLATION OF THE BELL-CHSH 
INEQUALITY 

The maximal value f max of the Bell-CHSH parameter 
dl·ll for an arbitrary mixed state was analyzed in Refs. fisl 
[Ï9J . For a pure state with concurrence C one has simply 
£ max = 2yT +C 2 H3. For a mixed state there is no one- 
to-one relation between C and £ max . Depending on the 
density matrix, £ max can take on vàlues between 2Cv / 2 
and 2vï + C 2 . The dependence of £ ma x on p involves 
the two largest eigenvalues of the real symmetric 3x3 
matrix R T R constructed from Rki = Tipat <8 07, where 
o\ = a x ,(T2 = Cy and 03 = a z . In terms of 71 and 72 , the 
elements Rm take the form 

(i + M 2 ) t t (1-M 2 ) t T 

Rki = V ^ ' ^ Tr 7l V fc 7i^ T + 1 ^ ' ; Ti-7 2 V fc72t 7 f T . 

(31) 

The matrix 72 has a polar decomposition 72 = Uy/ÇV 
where U and V are unitary matrices and £ is a diagonal 
matrix holding the eigenvalues of 7^72- The real positive 
£i's are determined by 

£i+6=Tr 72 7 2 , 2 v /^=|Tr 7 t7 2 |. (32) 
The matrix 7 i can be conveniently expressed as (see Ap- 



pendix [Bj 



7i = UQ^/lV, where Q = 




The parameters 01,02,03 are real numbers. The matrix Q 
is traceless due to the orthogonality of 71 and j 2 - The 
number c\ E (—1,1) on the diagonal is related to the 
inner product of 71 and 72 and takes the form 

Tr t 

ci = / 7l J 2 , with Tr7}7 2 = Tra z X^X. (34) 

Çl — Ç2 

The numbers 02,03 are determined by the norm and tilde 
inner product of 71 and satisfy the relations 

c 2 + c 2 c 3 = l, c 2 (ei + 6) + cÍ6 + c^i=Tr 7l t 7l . (35) 

We substitute 71 of Eq. and the polar decompo- 
sition of 72 in Eq. H31f) and parameterize 

3 3 
U^a k U = N ^ V(J l v ^ = Mki(T i' (36) 

t=l i=l 

in terms of two 3x3 orthogonal matrices N and M. The 
matrix R takes the form 

R = NR'M T , (37) 

where R! is given by Eq. (|31|l with substitutions R—>R', 
72 — > V? and 71 — > Qv 7 ?- With the hclp of Eqs. 
(|32I34I35|I . the eigenvalues it; of i? T i? can now be ex- 
pressed as (see Appendix lü|) 

U1 = 2A^( T+V/?JT ^)' (38) 



U2 = #( T -^)' (39) 



, |«| 4 |Tr 7l t 7il 
M 3 - 4 ^ 

where 



T = J\í 2 + MTtjIW 2 - 4(1 - H 4 ) (Tr 7l t 7 1 Tr7 2 t 72 - Tr 2 7Í7 2 ) » (41) 

V = 4|Tr 7l t 7l | 2 (Aí 2 - 4(1 - |a| 4 )Tr 7l t 7 1 Tr7 2 t 7 2 ) . (42) 

I 

We can relate the itj's to X^X and \a\ 2 using Eqs. (|14I23I25I26I34|I . The parameter £ max depends on the 
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FIG. 2: Parameter space of a beam splitter with [X^ X)u = 
1/2 spanned by \a\ 2 € (0,1) and |(X t X) H v| 2 € (0,1/4). 
All points correspond to a non-vanishing polarization- 
entanglement (C > 0) except the line segments \a\ 2 = and 
|(X^X)hv| 2 = 1/4 where entanglement vanishes (C = 0). 
Only in the shaded region, the Bell-CHSH parameter is able 
to detect entanglement (£ max > 2). The lines correspond to 
the functions /, g of Eqs. 1471481 respectively. 



two largest eigenvalues of R T R as 

£max = 2-v/ui +max(u 2 ,u 3 )· (43) 

Generically, the expression for £ max takes a complicated 
form where ordering of 112 and 113 depends on X^X and 

M 2 . 



V. DISCUSSION 

The objective of the discussion is to reveal the role 
played by the Mandel dip -2|a| 2 |(X t X) H v| 2 in the con- 
nection between C and £ max . 

We first consider the case |(X' i 'X)hv| 2 = 0. The con- 
currence of Eq. I|27|l reduces to 



C = 



2\a\ 2 Y{ l ^ y ^{X^X) lt {l~{X^X) ll ) 
(XtX)HH + (XtX)w - 2(XïX) Híi {XïX) 



vv 



(44) 

The maximal value of the Bell-CHSH parameter takes 
the form 



£ max = 2y/l +C 2 



(45) 



and C > implies £ max > 2. 

In the presence of a Mandel dip (\a\ 2 \(X^X) HV \ 2 > 0), 
the ability of £ to witness entanglement can disappear. 



We consider the special case (X*X)i 
currence of Eq. I|27|l reduces to 



|a| 2 (l-4|(XtX) 



HV 



l-4|a| 2 |(XtX) H v| 



1/2. The con- 



(46) 



To find £ max we have to consider the ordering of u-i and 
«3 which depends on |(XtJ s íT) H v| 2 and |a| 2 . The function 



/(H 2 ) 



2(1 + H 2 ) 



(47) 



divides parameter space in the region |(XtX)Hv| 2 < / 
where £ ma x = 2y/ui + U3 and the region |(XtX)Hv| 2 > / 
where £ max = 2-s/uï + uï- The equation fmaj; = 2 has a 
solution (?(|a!| 2 ) for |(XtX)Hv| 2 that lies in the region 
|(XtX)Hv| 2 < /■ The function g takes the form 

3(M 2 ) = \ (l - l«| 2 + l«| 4 ~ (1 ~ \a\ 2 )VT+W) 

(48) 

and breaks parameter space in two fundamcntal regions: 
a region \(X^X) RV \ 2 
|(XtX) HV | 2 > 9 whe 
regions in Fig. [21 The maximal value of the Bell-CHSH 
parameter is given by 



< g where £ max > 2 and a region 
e £ m ax < 2. We have drawn these 



£ max = 2C|arVl+M 4 
in the region |(XtX) H v| 2 < /■ 



(49) 



VI. CONCLUSIONS 

In summary, we have calculated the amount of 
polarization-entanglcmcnt (concurrence C) and its wit- 
ness (maximal value of the Bell-CHSH parameter £ ) in- 
duced by two-photon interference at a lossless beam split- 
ter. The ability of £ to witness entanglement (£ max > 2) 
depends on the Mandel dip -2|a| 2 |(X t X) H v| 2 · In the 
absence of a Mandel dip, C > implies £ max > 2 cf. Eq. 
(|45|l . whereas in its presence this is not necessarily true. 
In the latter case, as we have demonstrated in Sec. 
with (X'X)u — 1/2, the witnessing ability of £ depends 
on the individual contributions of temporal ( | et | 2 ) and 
polarization indistinguishability (|(XtX) HV | 2 ). 

Our results can be applied to interference of other 
kinds of partides, getting entangled in some 2® 2 Hilbert 
space and being "marked" by an additional degree of 
freedom. However, determining the indistinguishability 
parameter \a\ 2 requires careful analysis of the detection 
scheme. In case of fermions, the matrices 71 and 72 of 
Eq. (|TÏÏ|l are to be interchanged. Systems without a 
time-reversal symmetry are captured by the analysis, as 
we did not make use of the symmetry of the scattering 
matrix. 
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It is the purpose of this appendix to demonstrate that 
A = QVC where Q is a real traceless matrix of Eq. 133(1 
with Ci given by Eq. and c 2 ,c 3 satisfying Eq. (|35l) . 

The inner and tilde inner product of 71 and 72 take 
the form 

Tr 7+72 = 11 .4^, (B7) 



APPENDIX A: ARBITRARINESS OF 
TWO-PHOTON INPUT STATE 

The unitary scattering matrix has a polar decomposi- 
tion 



S 



K' 
L' 



K 
L 



> (Al) 



where K',L',K,L are 2x2 unitary matrices and T = 
diag(Tfj,Tv) is a matrix of transmission eigenvalues 
T H ,TV £ (0, 1). The outgoing state l^cmt) is related to 
the 4x4 matrix 



S 



CTin 




s 1 



(A2) 



cf. Eq. JSJ). By group decomposition K = K1K2 and 
L = L1L2, lí'out) is easily seen to correspond to K2<Ji n L 2 r 
scattered by S of Eq. (|A1|) with substitutions K — » K\ 
and L — > L\. 



APPENDIX B: JOINT SEMI-POLAR 
DECOMPOSITION 

The matrices 71 and 72 have a decomposition 

71 = UAV, 72 - Uy/ÇV, (Bl) 

where U ,V are unitary matrices and £ is a diagonal ma- 
trix of eigenvalues of 7^72- As we do not yet specify A, 
such a joint decomposition always exists. In our case, the 
matrices 71 and 72 have the special properties 



1*7i7i = -Tr 7 J,7 2 , Tr 7 {7 2 = 



(B2) 



TtjÍJíI = 2-JDet X^X Det(l — X^X), (B3) 



Tr 7Í7 X = Tr X^X - 2 Per X^X, (B4) 



Tr7Í7 2 = (DetUV)- 2 TrA^a y y/^cTy = 0. 



(B8) 



(Here we have used the identity Ua y U T = Det 2 Ua y , 
vàlid for any 2x2 unitary matrix U.) The conditions 
of Eqs. I|B7IB8(1 involve the diagonal elements of A as 
respectively 



Tr7Í7 2 = V^íi + V^m, 



(B9) 



Tr7Í7 2 = (DctUV)- 2 (y/frA^ + y^íi) = 0. 

(B10) 

It follows that Au = c í V6" an d A22 = — c*V^2 where Ci 
is given by 



Cl 



Tr 7i7 2 



£1-6 

The number ci is real since Tr7|7 2 = Tr <j z X^ X £ 



(Bll) 



The determinant of A is fixed by Tr'y\ > y 1 
implying 



Det.A = -y/&& 



-Tr7]7 2 



(B12) 



It follows that A12 = A'^e 1 ^ and A21 = -4 2 ie~ í0 with 
real A[2,A 2 i,(f>- The numbers -4' 12 ,-4 21 satisfy 



(B13) 



(B14) 



where Eq. i|B14(l comes from TiA^A = Trj\j 1 . The 
undetermined phase (f> can be taken out, 



A = 



e 2 







e" ! 2 



Au A'12 
A'21 A22 



e' 1 * 
e l í 



and absorbed in the unitary matrices U and V by the 
transformations 



Tr727 2 = TrX^X - 2~DetX f X, (B5) 



C7 



s *2 

e~ l í 



U. 



e~ l ~ 




U^U. (B16) 



Trjlj 2 = Tra z X^X. 



(Note that these transformations also hold for 72 since 
(B6) V? commutes with a diagonal matrix of phase factors.) 
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The matrix A is related to Q by A — Q\/£- It is now 
easily seen that the matrix Q is real and traceless and 
takes the form of Eq. QHHJI. with c\ given by Eq. 134(1 
and 02,03 satisfying Eq. 

As a last step we perform a consistency check to 
demonstrate that Eqs. (|B13IB14jl havc solutions for A' 12 
and -421- The Hermitian matrix X^X has an eigenvalue- 
eigenvector decomposition 



X^X = W^AW. 



(B17) 



^Í3 = ^(l + l«| 2 )ci(c 2 6 + c 3 Çi), (C2) 



R'22 = Jf(~ 1 + H 2 + í 1 + l«| 2 )( C ? + C 2 C 3)) 



(C3) 



In terms of the eigenvalues Aj € (0, 1) and the unitary 
matrix W, the inner product of 71 and 72 and the £j's 
take the form 

Tr7Í7 2 = Ai(|Wii| 2 - |W 12 | 2 ) + A 2 (\W 21 \ 2 - \W 22 \ 2 ), 

(B18) 



C^Atíl-Aa), 6=A 2 (1-A!). 



(B19) 



It follows that ci = cos 277, where we have set \Wi%\ = 
\W 22 \ = cos Tj and |Wi 2 | = |W 2 i| = sin??. Eqs. pi3IB14|) 
can be expressed as respectively 



A' 12 A' 21 = sin 2 2r? VAiA 2 (l- Ai)(l- A 2 ), (B20) 



4i = ^(l + N 2 )ci(c 2 + C3)v^6, (C4) 



^33 = ^ ((1 - M 2 ) + (1 + l«l 2 )c?) (6 + £2) (C5) 



The matrix iï' iï' has eigenvalues 



-^(l + l«| 2 )(c^2 + c^i). (C6) 



«1 



2AÍ 2 



T 



VT 2 - 42?) 



(C7) 



•A'u + ^'ai = siii a 2» 7 (Ai(l- A x ) + A 2 (1-A 2 )). (B21) 
Sincc 



"2 



2AÍ 2 



(C8) 



2VAiA 2 (l- Ai)(l- A 2 ) < Ai(l - Ai) + A 2 (l - A 2 ) 

(B22) 

a family of solutions exists. 



APPENDIX C: EIGENVALUES OF R T R 



The non-vanishing elements of R' are given by 



R' n = ^ (1 - |«| 2 - (1 + H 2 )(c 2 - C 2 C 3 )) y/Ç&, 



(Cl) 



u 3 — R'22 > 



(C9) 



where T,2? are the trace, determinant respectively of the 
2x2 real symmetric matrix 



R' n +R' 31 R'llR'l3 + R r 3iR's3 



( , j y iT " 3 

-R 11-R 13 + -R 31 -R 33 



i?' 2 , + i?' 2 



33 



(C10) 

By making use of Eqs. (|34I35|) us,T,T> can be simplified 
to yield the results of Eqs. I|4UI41I42II respectively. 
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